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Qh Abstract 

r^ Recent cosmological observations, including measurements of the CMB anisotropy 

"^ and the primordial helium abundance, indicate the existence of an extra radiation 

^ ^^ component in the Universe beyond the standard three neutrino species. In this 

paper we explore the possibility that the extra radiation has isocurvatrue fluctu- 
^^ ations. A general formalism to evaluate isocurvature perturbations in the extra 

P^ radiation is provided in the mixed inflaton-curvaton system, where the extra radia- 

\^ tion is produced by the decay of both scalar fields. We also derive constraints on the 

abundance of the extra radiation and the amount of its isocurvature perturbation. 
Current observational data favors the existence of an extra radiation component, 
but does not indicate its having isocurvature perturbation. These constraints are 
applied to some particle physics motivated models. If future observations detect 
^^ isocurvature perturbations in the extra radiation, it will give us a hint to the origin 

t> of the extra radiation. 
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1 Introduction 

It is known that most of the energy content of the present Universe consist of "dark" com- 
ponents: dark matter and dark energy. It is believed that the remaining components are 
well-known stuff: baryons, photons and neutrinos. These are ingredients of the standard 
cosmology supported by cosmological observations. 

However, there is no reason to exclude the existence of additional component other 
than listed above as long as its abundance is not large so that the cosmological evolution 
scenario is not much affected. Actually, the WMAP measurement of the cosmic microwave 
background (CMB) anisotropy combined with observations of the baryon acoustic oscil- 
lation (BAO) and the Hubble parameter (HO) constrains the effective number of neutrino 
species A^eff as Agfr = 4.34lQ;gg (68%C.L.) [1]. Including the Atacama Cosmology Tele- 
scope (ACT) data, the constraint is slightly improved as A'eff = 4.56 it 0.75 (68%C.L.) [2]. 
The recent results from South Pole Telescope (SPT), combined with WMAP, BAO and 
HO indicate A^eff = 3.86±0.42 (68%C.L.) [3]. Thus these observations suggest the presence 
of an extra radiation component beyond the standard three neutrino species at the 2a 
level ^. (See also Refs. [6, 7, 9, 8] for former analyses.) Hereafter we call non-interacting 
relativistic energy component other than neutrinos as "extra radiation" . 

The CMB measurement is sensitive to the expansion rate of the Universe, or the total 
energy density of the Universe, at around the recombination epoch. On the other hand, 
the big-bang nucleosynthesis (BBN) also gives a constraint on the expansion rate at the 
cosmic temperature of T ~ IMeV, where the weak interaction is frozen and the neutron- 
proton ratio is fixed. Thus the observation of the primordial helium abundance gives a 
constraint on the extra radiation component at T ~ IMeV. Recently, it was reported 
that the primordial helium abundance derived from the observations of the extragalactic 
HII regions indicates an excess at the 2a level compared with the theoretical expectation 
using the baryon number obtained from the WMAP data [10]. (See also Ref. [11].) 
Ref. [12] studied constraints on the effective number of neutrinos Agff and the mass of 
extra radiation component using CMB (including WMAP 7yr, ACBAR, BICEP and 
QUaD), SDSS and HO data, and similarly claimed the existence of an extra radiation 
and that its mass should be smaller than ~ 0.4 eV. If these observations are true, we 
need a theory beyond the standard model in which a new light particle species exists.^ 
Some cosmological scenarios and particle physics models are proposed in order to explain 
AA^eff ~ 1 [17, 18, 19, 20, 21]. 

The most literature treating the extra radiation so far (implicitly) assumes that the 
extra radiation has the adiabatic perturbation. However, the property of the fluctuation of 
the extra radiation depends on its production mechanism. For example, there is a scenario 
that the decay of a scalar field nonthermally produces an extra radiation component [17, 
20]. In this case the extra radiation has the fluctuation of the source scalar field, and 
if the scalar field is light during inflation, it obtains large-scale quantum fluctuations. 
Therefore, the extra radiation can have an independent perturbation from the adiabatic 



^ There are several studies which discuss how significant the deviation (A'off 7^ 3.04) is. For detailed 
discussion, we refer to Refs. [4, 5] 

^Discrepancy between the prediction and observation of the primordial helium abundance may be 
solved in the large lepton asymmetry scenario [13, 14, 15, 16]. We do not pursue this issue in this paper. 



one: it is an isocurvature perturbation. Investigating the isocurvature perturbation in the 
extra radiation may have a potential to distinguish the model of the extra radiation, if 
future observations further confirm AA'efr — 1. 

Therefore, in this paper we study the effect of the isocurvature perturbation in the 
extra radiation and derive constraints on them using the currently available datasets. 
The effect on the CMB anisotropy is similar to that induced by the neutrino isocurvaure 
perturbation. But we are mainly interested in the case where the energy density of the 
extra radiation itself causes significant change in the Hubble expansion rate. These effects 
of the extra radiation with isocurvature perturbation are phenomenologically taken into 
account by scanning both the effective number of neutrino species Nes and the neutrino 
isocurvature perturbation S^. 

First we develop formalism for calculating the primordial isocurvature perturbation 
in the extra radiation in Sec. 2. For concreteness we restrict ourselves to the two-scalar 
field case, which generalizes the curvaton model [22]. The extra radiation is assumed to 
be produced from decay of either or both of the scalar fields. We also point out that a 
large non-Gaussianity may exist in the isocurvature perturbation of the extra radiation. 
A similar study for the case of CDM/baryon isocurvature perturbations, including their 
non-Gaussianities, was done in Refs. [23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33]. In Sec. 3, 
we discuss the observational signatures of our model. In particular, we focus on the 
CMB angular power spectrum, which would be the best probe of our model for the time 
being. Then we derive observational constraints on such models with extra radiations 
and isocurvature perturbations using currently available datasets including WMAP, ACT, 
BAO and Hq measurement in Sec. 4. They are applied to some particle physics-motivated 
models in Sec. 5. We conclude in Sec. 6. 

We note that author in Ref. [34] also considered a specific model for isocurvature 
perturbations in extra radiation in the different context. 

2 Formalism 

We consider a cosmological scenario of two scalar fields, both of which obtain large-scale 
quantum fluctuations and may generate isocurvature perturbations as well as adiabatic 
ones. One is the infiaton field which initially dominates the energy density of the 
Universe during and soon after inflation. We also introduce another scalar field a. We 
do not specify which of them is the dominant source of the adiabatic perturbation, since 
it depends much on models. For convenience we call a the curvaton even if it may not 
dominantly generate the adiabatic perturbation. 

All the components of the Universe, photon, neutrino, CDM, baryon and extra radi- 
ation arise from both (p and a. In general there exist isocurvature perturbations among 
these components. But we particularly focus on the isocurvature perturbations between 
the extra radiation and the standard radiation. We denote the extra radiation component 
by X, photon by 7, and neutrino by z/. We also denote the Standard Model radiation 
which exists prior to the neutrino decoupling by r, and the X together with u by "DR" 
(dark radiation). We assume X is produced by the decay of and/or a and has no 
interaction with the other components throughout the whole history of the Universe after 



production. 

2.1 Nonlinear isocurvature perturbation 

We define ( as the curvature perturbation on the uniform- density slice. After the curvaton 
decay, it is the curvature perturbation on the shce where the total radiation energy density 
is spatially uniform.'^ According to the 6N formalism [35, 36], it is related to the local 
e-folding number N{tf,ti]x) as 

C{x) = N{tf,U;x)-ln^, (1) 

where a{t) is the background scale factor. Here the hypersurface at t = tf, in the radiation 
dominated era after the curvaton decays, is chosen to be the uniform density slice and 
the initial one at t = tj to be spatially flat slice. The curvature perturbation ( on large 
scales is conserved unless there are isocurvature perturbations, or non-adiabatic pressure. 
It is expanded by the scalar field fluctuations 6(j)i as 

C = N^M^ + ^N^.^.HMj + ■■■■ (2) 

In the two-field case, up to the second order, it is given by 

C = N^5(i) + N,5a + -N^^5(t)5(t) + N^M^^^ + -N^JaSa. (3) 

Similarly we define the curvature perturbation of the fluid i, Ci; as the curvature 

perturbation on the hypersurface where the energy density of i-th fluid is spatially uniform. 

Ci's are conserved on sufficiently large spatial scales as long as there no interactions or 

energy exchange between fluids [36]. If the i-th fluid has the equation of state pi = WiPi, 

(i and C are related by 

p,(x)=p,e3(i+-')K-C), (4) 

where Pi{x) is evaluated on the uniform density slice. Then we define the isocurvature 
perturbation between i-th and j-th fluids as [37] 

5,, ^3(0-0). (5) 

In this paper we are interested in the isocurvature perturbation in the extra radiation 
X. For convenience, we define the isocurvature perturbation in the dark radiation, which 
includes both X and neutrino, as 

^dr = 3(Cdr-C)- (6) 

It is formally expanded by the scalar field fluctuations 6(pi as 

5'dr = S^^S(f)i + -S^^^^S(j)iS(j)j + — (7) 

^For the reason discussed later, we distinguish the curvature perturbation at different epochs: (, ( 
and C- They correspond to the curvature perturbation after the curvaton decay, neutrino freezeout and 



p± 



annihilation, respectively. It is C that is directly related to cosmological observations. 



In the two-field case, up to the second order, it is given by 

5'dr = S^5(j) + S^6a + -S^^6(f)64> + S^JcpSa + -S^^6a6a. (8) 

The power spectra of the curvature/isocurvature perturbations, and their correlation 
are given by 

(C(^i)5dr(^2)) = {27,f5{k, + 4)P^Sbr,(^i), (9) 

where 

Pa{k) = [Nl + N'„]Ps^{k), 

PcsnM = [N<pS^ + N,S^]Ps^{k), (10) 

Ps^^s^Ak) = [Sl + Sl]Psm. 
where we have neglected higher order terms, and the power spectrum of S(j)i is defined as 

{S(l)i{h)S(l)j{k2)) = {27Tf6{h + k2)P54k,)6,,, (11) 

Here ifjnf is the Hubble parameter during infiation, rig is the scalar spectral index^ and 
fco is the pivot scale chosen as k^ = 0.002Mpc~^. 

The correlation parameter between the curvature and isocurvature perturbations is 
defined by 

^ P(SuB.if^o) ^ N^S^ + N^S^ 

~ VPcdko)PsunSuniko) ^[Nl + Nl){Sl + Sl) 

The uncorrelated isocurvature perturbation corresponds to 7 = 0, and totally (anti- 
) correlated one to 7 = (— )1. The effect of isocurvature perturbations on the CMB 
anisotropy depends on its magnitude as well as the correlation parameter. Thus what we 
need is to express A^<^, S^, and so on, in terms of model parameters. 

For later use, we also define the dimensionless power spectrum as VAB{k) = {k^/27i'^)PAB{k) 
where the subscripts A and B are either ( or S'dr. Using the quantities given above, they 
are expressed as 

ncW^K + iv,T(t^)^(^)"". 

T<s^M = mS* + iV.S,l (1^)' (i)"* ' , (14) 

^w„w^K+sji(|f)^(i)"^". 

In the following we formulate a method for calculating the extra radiation isocurvature 
perturbation in a typical cosmological setup. 



^ The scalar spectral indices for </) and a do not coincide in general. In the following we assume they 
are the same just for simplicity. 



epoch 


component 


energy transfer 


r^<H 


0, a 


-^ XW + rW 


T^<H <T^ 


XW, rW, a 


a -^ X('^) + r^'^) 


T,<H <T^ 


X,r 


r — )■ Z/ + Te 


Te± <H < r^ 


X, u, re (DR = X + u) 


e=^ -^ 7 


H <T^± 


X, J/, 7 (DR = X + u) 





Table 1: Energy components of the Universe at each epoch except for CDM and baryon. 
Here T^lTcr) is the decay rate of the inflaton (curvaton), and T^, denotes the neutrino 
interaction rate at the neutrino freezeout. Fei denotes the Hubble parameter at the e^ 
annihilation. Here Tg denotes the plasma consisting of 7 and e^. 



2.2 Extra radiation and isocurvature perturbation 

Now let us move to our concrete setup. Here we consider as the inflaton field. We want 
to know the final dark radiation ("DR") isocurvature perturbation, and hence we must 
connect it to the primordial perturbations of and a through some transition points 
where constituents of the cosmological fluids change. Thus we evaluate the curvature 
perturbations step-by-step in the following. The evolution of the fluids in our model is 
summarized in Table. 1. We assume that the a decays well before BBN begins. Con- 
sidering that the decay of a after BBN might produce significant energy in the form of 
radiations and hadrons and upset the standard BBN, it is a reasonable assumption. In 
the most part of this paper we take this assumption.^ 

In general, the CDM and baryon can have isocurvature perturbations depending on 
their origins. However, including them makes the analysis too complicated. Thus we focus 
on the case where the extra radiation has an isocurvature perturbation and CDM/baryon 
do not. For example, if CDM/baryon are created from thermal bath after the curvaton 
decays, they only have adiabatic perturbations and our assumption is justified. 



2.2.1 At the curvaton decay 

Let us take the uniform density slice at the curvaton decay : H{x) = T„, where F^- is 
the decay rate of a. We make use of the sudden decay approximation in the following 
analysis [38]. On this slice, we have following relations 



S'i') 



py{x) + {l- r„)p^{x) = px{x), 

p^^\x) + r„p„{x) = Pr{x), 

pi^\x) + pf{x) + p^{x) = ptot(f)(= Ptot), 



(15) 
(16) 
(17) 



where r„ is the branching fraction of a into particles other than the extra radiation X. 
Superscript (0) means that the corresponding component comes from the inflaton decay. 



^ The BBN constraint can be avoided if the a decays dominantly to X particles even if its decay 
occurs after BBN. For completeness, we will also discuss such a case in Appendix. 



The curvature perturbation of each component is related to its background value as 









Pre'^'-'\ 



Px{x} 

Pr{x) 
Pa{x) 

From these equations we obtain 

i?«')e4(C.-C) + ^^ 
i?We4(C0-C) + 4-)e3(C.-C) = R^e^iCr-0 ^ 

(1 - i?,)e^(^^-^) + R^e'^^''-^^ = 1, 



,{^)g3{C.-C) 



i?;.e^(^--«, 



where we have defined 
and 



Ru = Pa/ptot at o- decay, 



(18) 
(19) 
(20) 
(21) 
(22) 

(23) 
(24) 
(25) 

(26) 



R 



Px /Ptot 



X 
X 



Px Vp 



4^^ - p1^Vp-- 



tot 
tot 
'tot 



(l-r^)(l-i?.), 
(1 - r,)R^, 

r^(l --Ra), 

r R 



(27) 



where r^ is the branching fraction of into particles other than the extra radiation X,^ 

Pr/Ptot- Notice that Rx + Rr = 1- 



M 



(^) 



Ri^^ + Ri^^ 



and Rx = RT + Rx = Px/Ptou R 
All of these quantities are evaluated at the curvaton decay, H = T^. As is already stated, 
these quantities are not constants in time because of the changes in the relativistic degrees 
of freedom g^. Solving these equations up to the second order in ^^ and (a, we obtain 



c 

Cx 



1 



+ i?(C - C<^) + ^R{i -R){S + R){C^ - C 



4R 



X 



{AR 



W 



X 



(3 + R)R 



X 



8Rx 



{1-R)rP)Q + {3 + R)rP(:, 
[1 + R){3- R)Rf + (1 - R)RR 



X 



(a - U' 



(28) 



C, = -^ [(4i?W + (1 - i?)i?('^))a + (3 + R)R^:K.] 

+ ^^ ^«^2^^''^ [(1 + ^)(3 - R)R^J'^ + (1 - R)RR^:^ iU - Ca) 



8RI 



^Precisely speaking, 1 — r^ is not the branching ratio of (j) into X (B^^x)- It is given by 1 — r^ 



B. 



^^xe/{l - B^^xil - e)) where e = [g,iH = T^)/g,{H 



11/3 



with r^ denoting the inflaton decay 
width, taking into account the change of relativistic degrees of freedom between the inflaton decay and 
curvaton decay. 



where we have defined 

Note that Cx and Cr are conserved quantities at superhorizon scales [37]. In general, the 
curvature perturbation ( is not conserved in the presence of non-adiabatic pressure. In 
the present case, both X and r behave as radiation, and hence one may consider that 
( is conserved on sufficiently large scales. More precisely, however, the radiation energy 
density pr does not scale as a(t)~'^ because of the changes in the relativistic degrees of 
freedom g^. Therefore, ( is not conserved even after the curvaton decay. The changes in 
( at the changes in g^: depends on the fraction in which the X dominates the Universe. 

2.2.2 At the neutrino freezeout 

Next let us take the uniform density slice at the freezeout of neutrinos H = T^ {T ^ 
IMeV). On this slice we have the following relations 

(1 - Cu)pr{x) = p-y{x), (30) 

CuPr{x) = Pu{x), (31) 

Prix) + pxix) = ptot(f)(= ptot), (32) 

where Ci, is the energy fraction of neutrinos. For the standard three neutrino species, it 
is given by 

c,y = — = — at neutrino freezeout. (33) 

Here the symbol 7 should be interpreted as a thermal plasma consisting of photons elec- 
trons and positrons. The curvature perturbation of each component is related to its 
background value as 

p,{x) = p,e'^<'^-^\ (34) 

p^(f)=p^e^(C.-C), (35) 

where ( denotes the curvature perturbation at this epoch and it is different from (^ in 
Eq. (28). From Eqs. (30) and (31) we find 

C7 ~ Cr; Qi^ = Cr; (3DJ 

which holds nonlinearly, as is expected. We also define dark radiation (DR) as the sum 
of extra radiation X and neutrino. It satisfies 

pxix) + Cyprix) = pxix) + puix) = Pdr(^), (37) 

on the uniform density slice at the neutrino freezeout. The curvature perturbation of DR, 
Cdr, is given by 

PDRix) = PDRe"(^°«-^^ (38) 



From Eq. (37), we find at first order 



<,DR — SX + 'sv 

Pdr Pdr 



-Kdr -Kdr 



(39) 



where we have defined i?DR = -Rdr + -^dr = PDR/ptot|ii-=r,, with 



-TLriR — ^X ' ^i^-^r ' 



i? 



>(0) 
^-DR 

DR 



iXv ~r Ci/SX:^ 



(40) 



^-x 



Here quantities with a tilde means that they are evaluated at the neutrino freezeout. 
Relations between quantities with and without a tilde are given by 



R 



R 



X 



X 



Dip g*{H=T„) 



1/3-' 



Rr 



p / g4H=Ta) 
^^ \a,(H=T,) 



1/3 



p I p g*{H=V„) 



1/3- 



(41) 



This is approximated as 



R 



X 



R 



X 






1/3 



J. tjip J. Lr* 



(42) 



if px{H = Tu) ^ Pr{H = T^). Here g*{H = F^) = 10.75 is the relativistic degrees of 
freedom at the neutrino freezeout. Using these quantities, we find that ( is given by 



+ 



3 + R ( RrRi"^ 



4 \ R, 

{C-Q + 



+ 



RxRx 



R 



X 



1 + R){3-R) 3 + RlRrR^^ RxRx\\,. ..2 

(C<7 - Q^) 



Rr 



+ 



R 



X 



(43) 



It can be checked that it coincides with C^ in Eq. (28) for Rr = Rr and Rx = Rx- 
Therefore, as is expected, the relative difference between the evolutions of px and pr 
results in the non-conservation of the curvature perturbation. 

Up to the second order in ^^ and („, the DR isocurvature perturbation is given by 



5'dr = 3(Cdr — C) 



3 -\- R RrRx I R 



{^) 



R 



X 



^ -RdR V ^r 



R 



X 



(1 - c.) 



5M/ 



-(c - a) 



A)/ 



^ {R-3){R+l) ^ 3J^ / RrRr\c, + Rj,^) ^ RxR'x\l + Rt,k) 



2-Rdr 



Rr 



R 



X 



(C - C^)^ 

(44) 



2.2.3 After the e^ annihilation 

Moreover, shortly after the neutrino freezeout, e^ annihilation takes place. This also 
slightly affects the curvature perturbation, since the radiation energy increases while other 
components {X and z/) are unaffected. By taking the uniform density slice well after the 
e^ annihilation, we find 



Py{x) + Pu{x) + px{x) = Ptot{x){= ptot)- 

In a similar way, we find the curvature perturbation ( at this epoch as 



(45) 



c 



+ 



3 + R { RrRi"^ , RxRx^ 



4 \ R, 

(C-Q + 



+ 



R 



X 



1 + R){3 — R) 3 + R I RrRr RxRx \ { //■ /■ \2 



2 2 \ Rj. Rx 

where quantities with a hat are defied at this epoch. They are given by 

Rx f, (1 + S)Rr 



(46) 



R 



X 



where 



'l + 5)Rr + Rx 

1/3 



Rr 



;i + 6)Rr + Rx 



1 + 6 



[1 - Cu) + C 



U I 1^ '^v 



1.205. 



(47) 



(48) 



Here the "r" should be interpreted as the sum of 7 and v. The result (46) has the 
same form as (43) except for quantities with a tilde are replaced with those with a hat. 
Since Cdr is conserved through the e^ annihilation process, we can also evaluate the DR 
isocurvature perturbation as 



'S'dr = 3(Cdr — C) 



3 -\- R RrRx I R 



{^) 



R 



X 



^ -RdR V ^r 



R 



X 



'1 - c„ 



5M/ 



-(C - Q 



»/ 



^ {R-3){R+l) ^ 3J^ / RrRi^>{c, + R^n) ^ RxRFi^ + Rdr) 



2-Rdr 



Rr 



R 



X 



(49) 



where we have defined -Rdr = PoR/ptot evaluated after the e^ annihilation, which is given 

by 



and 



P^ 



P1/ + P7 



-Rdr — Rx + CjyRr, 

0.405 after e^ annihilation. 



(50) 
(51) 



Notice that £^{1 + S) = Cy. 
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2.2.4 Connecting to primordial perturbations 

Finally we relate the curvature perturbation of and a to their quantum fluctuations, 
6(j) and 6a. Let us take the uniform density slice when the curvaton a begins to oscillate. 
Assuming that the curvaton energy density at this epoch is negligible and the total energy 
density is dominated by the infiaton or the radiation coming from the infiaton decay, the 
curvaton oscillation begins uniformly in space on this slice. Just after that, the curvaton 
behaves as matter.^ Thus we have 



p^[x) = p„e 



3(Ca-C<^) 



(52) 



From this we obtain 



c. 



1 /M._15pl 

3 V p. 2 p^ 



1 f25a 

3 V CT,; 



5a^ 



(53) 



where cxj denotes the initial amplitude of the curvaton during inflation. On the other 
hand, C,^ is given by the standard formula 



1 ^.. 



1 



M^V, 



2Ml 



VI 



m' 



(54) 



where V is the infiaton potential and V^iV^fj,) are its first (second) derivative with respect 
to 0. 

Now we have obtained all ingredients for expressing curvature/isocurvature perturba- 
tions in terms of model parameters. From Eqs. (46), (54) and (53), we obtain 






1 V 



X 




3 + AR-2R'-2{3 + R) 



rifR 






/ J 



Rr 



R 



(55) 



^ We assume the curvaton has a simple quadratic potential for simplicity. Otherwise, analytic estima- 
tion is difficult particularly for the analysis of non-Gaussian fluctuation [39, 40]. 
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From Eqs. (49) and (53), we obtain 



S. 



S. 



3 + XL R^Rx I R 



H 



R 



X 



2(Ti R 



DR 



Rr 



R 



x 



3 + XL Rj-Rx I Rr 



(<^) 



R 



x 



2^' i?DR V Rr 

2R^ - 4i? - 3 

Sa, = Srhrh = 0. 



R 



X 



(1-c.) 
3 



(56) 



?M^ 



2( 3 + i?) / RrRl-'"{c, + RBR) _^ /^xi^^^(l + ^DR) ^ 

^DR \ Rr Rx J 



For later use, we summarize quantities relevant for calculating the CMB anisotropy. 
We define the effective number of neutrino species, Nes, by 



Pdr = P,. + Px = A/eff- I — 1 p^. 



(57) 



where quantities are evaluated after the e^ annihilation. The extra effective number of 
neutrino species, ANgg, is then given by 



AN, 



px SRx 3Rx 43 i? 



X 



eS 



Pv 



C-^ijjrLj' L^y-Lhy jrij'p 



(58) 



While our primary interest in this paper resides in the Gaussian perturbations and hence 
the primordial power spectrum, the primordial bispectrum generated from non-Gaussian 
perturbations will be briefly discussed in the next subsection. The isocurvature pertur- 
bation in the dark radiation is then given by 
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The curvature perturbation is given by 

C = N^6(f) + NM = 



+ 



3 + R ( RrRi"^ 
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Rr 
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R 



X 



CTi 



(60) 



The correlation parameter (13) in this case is given by 

7 = sign(S'<,^ 



m + m 



(61) 



2.3 Non-Gaussianity in isocurvature perturbations of the extra 
radiation 

In this paper we do not go into details of analysis on non-Gaussianities, but it may be 
worth mentioning them briefly here since, to the best of our knowledge, no literature 
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has discussed non-Gaussianities in the extra radiation isocurvature perturbation or the 
neutrino isocurvture perturbation. For this purpose, we have expressed all quantities up 
to the second order in 50 and 5a. 

Non-Gaussianities of the cosmological perturbations are characterized by their bispec- 
tra. In the present case we have isocurvature perturbations in the extra radiation as well 
as the adiabatic perturbation, and hence non-Gaussianities may appear in various combi- 
nations. We define the bispectrum of the curvature/isocurvature perturbations through 
their three point correlation functions as 



(62) 



UkiK{k2)ah)) = {27rfS{h + h + h)B^cdki, h, h), 
{ah)ak2)SMh)) + (cyclic) = {27T)'6{h + h + h)Bas{ki, k^, fca), 
(C(^i)^dr(^2)^dr(4)) + (cyclic) = {27Tf6{h + k2 + h)Bcss{ki, k2, fcs), 
{SDR{ki)SnR{k2)SBR{h)) = {2Tif6{h + k2 + h)Bsss{ki, k2, fca), 

where 

Bccdku k2, ks) = [N^N^^ + N'^N^^] [P5^{h)Ps^{k2) + (2 perms.)] , 

%5(fci, k2. h) = [NlS^^ + NlS,„ + 2N^S^N^^ + 2N,S„N,„] [Psdki)Psdk2) + (2 perms.)] 

Bcssiku k2, ks) = [SlN^^ + SJN^^ + 2N^S^S^^ + 2N^S^S^^] [Ps4>{ki)P54>{k2) + (2 perms.)] , 

Bsssiki, k2, ks) = [SlS^^ + SJS^^] [P&dki)Ps<t>{k2) + (2 perms.)] . 

(63) 

Here we have neglected contributions from purely non-Gaussian parts, e.g., terms which 
are proportional to A^^^, N'^^S^^, etc. The first one, -Biji^i^, is same as that usually studied 
in the mixed infiaton-curvaton system [41]. Others are bispectra that arise only when 
the extra radiation has isocurvature perturbations. Also one should notice that the first 
two terms and the remaining two terms in -B^^s and Bc^ss have different effects on the 
CMB anisotropy. Non-Gaussian imprints of the extra radiation components on the CMB 
anisotropy may also be important for constraining or confirming models of extra radiation. 
We leave these issues for future work. 

3 Observational signature of isocurvature perturba- 
tions in the extra radiation 

Isocurvature perturbations in the extra radiation generated in the early Universe affect 
the late-time structure formation and hence the CMB anisotropy. Since extra radiations 
and active massless neutrinos are cosmo logically equivalent, the initial condition for the 
structure formation with S'dr 7^ can be identified as the neutrino isocurvature pertur- 
bation mode [42]. In this section, the neutrino isocurvature density (NID) perturbation 
collectively means the isocurvature perturbation of the dark radiation (DR), which is sum 
of the active neutrinos and the extra radiation particle X. Therefore, all the following 
analyses are applied to the case of the standard neutrino isocurvature perturbation if 
there are no extra radiation particles and the dark radiation consists only of the active 
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neutrinos: ^dr = 5'y(= 3(Ci/ — C)) and i?DR = Ru{= Pu/Ptot)- Our formulation includes 
more general case where the extra radiation X significantly contributes to the relativistic 
energy density and the isocurvature perturbation. In addition, throughout this section, 
we assume a vanilla fiat ACDM model as the unperturbed Universe, and all the associ- 
ated cosmological parameters are fixed to the WMAP 7-year mean value [1], except the 
effective number of neutrino species, which can deviate from the standard value due to 
the existence of the extra radiation, i.e. N^s 7^ 3.04. 

So far, several authors have investigated the initial perturbation evolution (i.e. at 
superhorizon scales deep in the radiation dominated (RD) epoch) for the NID mode [42, 
43]. Perturbation evolution at subhorizon scales or in later epoch can be calculated 
numerically by adopting Boltzmann codes (e.g. CAMB [44]). This allows us to constrain 
the amplitude of the NID mode from current cosmological observations. For the time 
being, the best probe of the NID mode is the CMB temperature anisotropy. Thus we 
concentrate on the CMB signals of the NID mode in this section. 

Although the adiabatic (AD) and the CDM isocurvature (CI)*^ modes have been inves- 
tigated in many literatures (e.g. Ref. [47]) and well understood, there are few papers which 
offer physical understanding of the perturbation evolution beyond the initial behavior for 
the NID mode (See also Ref. [48]). In the following, we investigate the semi-analytical 
description of the perturbation evolution for the NID mode. 

3.1 Perturbation equations 

We firstly summarize the evolution equations of the cosmological perturbation at a linear 
level. Scalar perturbations in the fiat Universe are to be discussed, which can be expanded 
in terms of the eigen function Q{k, x) = ex]i[ik-x]. We consider a single Fourier mode with 
a wave number k. In addition, we adopt the conformal Newtonian gauge, where physics 
can be intuitively understood in various cases. In this gauge, the perturbed metric is 
given by 

ds'^ = a{r]f [-(1 + 2'^Q)dri^ + (5^(1 + 2^Q)dx'dx^] , (64) 

where 77 is the conformal time, airf) is the scale factor, and \1/ and $ are the gravitational 
potentials. The perturbed energy-momentum tensor of a fiuid is given by 

T° = -(l + 5Q)p, (65) 

Tf = VQ,{p + p), (66) 

TJ = -VQ^{p + p), (67) 

T; = -{S]{l + 7iLQ) + UQi)p, (68) 

where 6 = 6p/p, V, tcl = Sp/p and 11 are the density, bulk velocity, pressure and 
anisotropic stress perturbations, respectively. Note that Qi = —ikiQ and Qij = —{kikj — 
Sij/3)Q represent (traceless parts of) the spatial derivatives of Q. Our notation is basically 
same as that adopted in Ref. [47]. 



*We omit the baryon isocurvature mode since its CMB power spectrum completely degenerates with 
one for CI (See e.g. Refs. [45, 46]). 
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The Einstein equation gives a set of perturbation equations for the metric, 

k^^ + m{^-H^) = 47rGa2pt„Aot, (69) 

k{^-H^) = _4vrG'a2(ptot+Ptot)Kot, (70) 

$ + 7^(2<i.-^)-(2-H + -H2)^ + |-($ + ^) = 47rGa2ptot VTL tot, (71) 

A;2($ + ^) = -8nGa%,,U,,t (72) 

where "H = a/a is the conformal hubble expansion rate and the subscript "tot" represents 
the total fluid which consists of CDM, baryons, photons and neutrinos. Here and hereafter 
the dot represents a derivative with respective to rj. 

On the other hand, the conservation law of energy-momentum gives the perturbation 
equations of the fluids. For CDM they are given by 

S^ = -kVc - 3$, Vc = -HVc + A;^. (73) 

As for baryons, we obtain 

K = -kV,-m, H = -^H + A;^ + ^^(V;-H), (74) 

where / = 3ph/4:p^ roughly denotes the ratio of the baryon and photon energy density 
and the cross section of the Thomson scattering and the number density of electrons are 
denoted by cxr and n^, respectively. We also obtain perturbation equations for photons 

5^ = -^kV.,-4:^, V^ = k{5^ + <iI-^n^)-aaTneiV^-Vb), (75) 

o u 

and DR 

<^DR = -^fcVDR-4<i>, \/DR = fc(5DR + ^- JHdr). (76) 

3 o 

Note that here we have omitted evolution equations for the anisotropic stress and higher 
order multipole moments of photons and DR, which are not very important for our dis- 
cussion. (However, these are included in numerical calculations.) 

3.2 Initial conditions 

By solving equations (69)-(76) at superhorizon scales deep in RD epoch, we obtain several 
independent initial solutions for the structure formation, which can be characterized by 
the gauge-invariant curvature perturbations ( and isocurvature perturbations Sij ^. While 
our primary interest resides in the NID mode, other initial modes (i.e. AD and CI) are 
also discussed in a parallel manner for reference. 

Now let us summarize the initial conditions for the AD, CI and NID modes. In 
Ref. [42], they are derived in the synchronous gauge, and we can obtain the initial condi- 
tions in the Newtonian gauge by performing a gauge transformation. Since we are mostly 



^ In general, there can be a neutrino velocity (NIV) mode, which cannot be characterized by C nor 
Sij [42]. We do not consider the NIV mode in this paper. 
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interested in the CMB signals of these modes, here we restrict ourselves to the initial 
solutions for only the perturbations in the photon fluid and the metric. For later conve- 
nience, we define fceq = flmHo/y/Tl^., which roughly corresponds to the horizon scale at 
the matter-radiation equality. 

The AD mode is characterized with nonzero initial (. Up to 0{r]), the initial condition 
for the AD mode is given as 

^,/C(o) ^ 'r + Hi^ + i^R.n)^ 

15 + 4i?DR 2(15 + 4/2DR)(15 + 2i?DR) 
V,/m = - ^^ ]^ kv, (78) 

15 -h 4i?DR 

»/C(0) = ^_ + 2,i3-SR.,) 

15 + 4flDR 8(15 + 4flDE)(15 + 2flDR) 

*/C(o) = H±i|25 ^I|±i5;?2i0_,,,„, (80) 

15 + 4flDR 8(15 + 4flDR)(15 + 2flDR) 

where -Rdr = - 'T^ is the fraction of DR (= neutrinos + X) in the radiation energy 
density, which can be rewritten in terms of A'^efr defined in Eq. (57) as 

We note that the initial density perturbation of photons 6^ and the metric perturbations 
\E', $ are all nonzero and comparable in the amplitude. As is well known, this fact results 
in the cosine-like acoustic oscillation of the photon-baryon fluid, which will be discussed 
in more detail in the next subsection. 

The CI mode is characterized by a nonzero initial CDM isocurvature perturbation 
Sc = 3(Cc — C) = 3(5c — ^S-y), where (c is the curvature perturbation on the uniform 
density slice of CDM. In the same way as the AD mode, the initial condition for the CI 
mode can be obtained as 

A /crm 4(15 + 4Rdr), ..^x 

6,/SM = -2(15 + 2^^^)^-., (82) 

V,/S,iO) = 0{if), (83) 

where Re is the fraction of CDM in the energy density of matters, i.e. Re = Pc/{Pc + Pb)- 
In the CI mode, 5^, "^ and $ all vanish initially, which results in the sine-like acoustic 
oscillation of the photon-baryon fluid. 
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The NID mode is characterized by a non-vanishing initial DR isocurvature perturba- 
tion S'dr. The initial condition for the NID mode is given as 

i?,(15 + 4/2DR) (2/?DR + 15)(4i?DR+15) '" 

V,/SMO) = ^ , °^ -A:r7, (87) 

15 + 4i?DR 4(4i?DR + 15)(2i?DR+15) 

$/^dr(0) = -^^ RM2R^K - 75) 

4i?DR+15 4(4i?DR + 15)(2i?DR+15) '" 

where i?^ = 1 — -Rdr- This NID initial condition is distinct from both the AD and 
CI modes in several respects. First of all, photons have non-vanishing initial density 
perturbations 5^(0) 7^ in the NID mode, which is in contrast to the CI mode, where Sy 
vanishes initially. On the other hand, while 6^ does not vanish initially for the AD mode 
also, the gravitational potentials for the NID mode are initially much smaller than those 
for the AD mode. This is because the density perturbations in photons and DR cancel 
out each other for the NID mode, and the density perturbations in the total radiations 
vanish. However, due to the presence of the anisotropic stress in DR, initial \E' and $ 
do not vanish exactly. In a nutshell, the NID initial condition can be characterized by 
non- vanishing 6^ and small $ and \E'. 

3.3 Acoustic oscillation 

Now we discuss how evolution of the photon perturbation in the NID mode differs from 
those in the AD and CI modes. Before recombination, as perturbations enter the horizon 
A; ~ "H, the photon fluid undergoes the acoustic oscillation. Tight-coupling approximation 
is a good approximation above the diffusion scale of photons. In the limit of tight- 
coupling, the photons and baryons can be treated as a single fluid with sound velocity 



1/ V 3(1 -|- /) and obey equations of motion given by 



6^ = -^A;V;-4<i>, (90) 

% = -^V, + ^kc% + k<iI, (91) 



which are combined into a single second order differential equation 

/ i , h.2„2s 



Sj + Y^^f + ^^''^7 = 4^> (92) 



where 



F = -^-^^--^. (93) 
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Figure 1: Perturbation evolution of the photon fluid and the metric. We showed the 
metric perturbations \E' (solid red) and $ (dashed green), the photon effective temperature 
perturbation S^/A + \E' (dotted blue) and the photon velocity perturbation V^/A/3(dot- 
dashed magenta). Positive and negative values are represented by thick and thin lines, 
respectively. From the top to bottom, each row corresponds to AD, CI and NID modes, 
in order. From the left to right, each column corresponds to three different scales k = 1, 
0.1, 0.01 Mpc-i in order. 

The left hand side of Eq. (92) describes the acoustic oscillation of the photon-baryon fluid 
and 4F in the right hand side can be regarded as an external force. 

Fig. 1 shows perturbation evolution of the photon fluid and the metric for the AD, CI 
and NID modes at three different scales [k = 1, 0.1 and 0.01 Mpc""^), which are calculated 
by using a modified version of CAMB. In the figure, we plotted the metric perturbations 
\E' and $, the photon effective temperature perturbation (5^/4 + ^ and the photon velocity 
perturbation V.y/\/3. As can be see from the figure, when perturbations enter the horizon, 
the photon fluids starts the acoustic oscillations, and at subhorizon scales the effective 
temperature S^/4 + ^ keeps oscillating with an almost constant amplitude as long as it 
damps inside the diffusion scale. On the other hand, the photon velocity perturbation 
V^/^/S oscillates with the same amplitude but a phase shifted by 7r/2 compared with 
the effective temperature perturbation. We also note that gravitational potentials decay 
during the horizon crossing in RD. 

Now we discuss the differences in the photon acoustic oscillation among the three 
different initial modes. The most distinct one that is found in Fig. 1 would be the difference 
in the oscillation phase for CI mode against other two modes. This originates from the fact 
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that in CI mode oscillation starts with vanishing 6^ at rj = 0, while oscillations of the other 
two modes start with no n- vanishing 6^. On the other hand, the differences between AD 
and NID modes may be less noticeable as their oscillation phases look similar. The most 
significant difference can be found in the change in the oscillation amplitudes between the 
horizon crossing. The amplitude of the acoustic oscillation gets boosted by the decay of 
gravitational potentials in the AD mode, which is less effective in the NID mode due to 
the initial smallness of the gravitational potentials. 

In order to verify the above discussion, we adopt the semi-analytic solution of the 
acoustic oscillation derived in Ref. [47]. Under the adiabatic approximation, the WKB 
solution of Eq. (92) can be obtained as 



[l + fiv)fU,iv) = S,iO)coskrM+'^ 



k 



KiO) + \f {0)5,(0) 



sin kr sit]) 



where 



+ ^ f dv' [1 + fiv')f'sm [kr.ir]) - kr^)] Fiv'), (94) 

rs{r]) = / dr]'cs{r]') (95) 

is the sound horizon of the photon-baryon fluid. In Eq. (94), the first and second terms 
correspond to the free damped oscillation which reflect the initial condition of photon 
density and velocity perturbations. On the other hand, the third term corresponds to the 
forced oscillation driven by the gravitational potentials. In the following discussion, the 
first, second, and third terms are called the cosine, sine and forced terms, respectively. 

In Fig. 2, the time evolution of the photon density perturbation are plotted for the AD, 
CI and NID modes at /c = O.lMpc"^. In the figure, we have shown exact numerical solu- 
tions for 5, as well as the WKB solutions Eq. (94) obtained by adopting the numerically 
evaluated source term F. We can see that the WKB solution is in excellent agreement 
with the exact one, regardless of initial conditions. This allows us to understand correctly 
the physical origin that drives the evolution of photon perturbations by comparing the 
different terms in Eq. (94). For this purpose, we have also plotted the WKB cosine, sine 
and forced terms separately in the figure. 

For the AD mode, we see that the forced term dominates the acoustic oscillation and 
imprint of the initial cosine-like oscillation induced by (5^(0) is almost invisible. However, 
the external force also induces an acoustic oscillation with cosine-like phase, as the gravi- 
tational potentials also oscillate in a similar manner around the horizon crossing. Totally, 
the acoustic oscillation for AD mode can be regarded as a kind of forced oscillation with 
a cosine-type phase. 

In the case of CI mode, the source term is also dominant as can be seen in Fig. 2. 
Since the gravitational potentials grow proportionally to r] at superhorizon scales as the 
density perturbations in CDM generates the curvature perturbation, they induce a sine- 
type acoustic oscillation rather than a cosine-type one. The acoustic oscillation in the CI 
mode can be regarded as a forced oscillation with a sine-type phase. 

On the other hand, in the case of NID mode, the cosine term dominates the acoustic 
oscillation and the forced term is subdominant. This is in contrast to the previous two 
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Figure 2: Time evolutions of the photon density perturbation S.y for AD (top), CI (middle) 
and NID (bottom) modes at k = 0.1 Mpc~^. Plotted are the exact solution (solid red) the 
total WKB solution (dashed green), and three terms consisting the total WKB solution 
including the cosine term (dotted blue), sine term (dot-dahsed magenta) and term from 
the external force (dot-dot-dahsed light blue). Positive and negative values are represented 
by thick and thin lines, respectively. 
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modes. This result comes from the fact that the gravitational potentials are initially 
small and grow little before the horizon crossing in NID mode. Therefore, the acoustic 
oscillation for the NID mode can be regarded as free oscillation with a cosine-type phase. 
Smallness of the gravitational potentials also affect the perturbation evolution of pho- 
tons after the recombination, when photons become coUisionless and free-stream in the 
potentials which are dominantly formed by the CDM density perturbations. While grav- 
itational potentials is constant deep in the matter dominated (MD) epoch, it is not fully 
MD just after the recombination, where the potentials decay during the horizon crossing. 
Effective temperature of free-streaming photons are boosted (or suppressed) due to the 
decay of the potentials $ — \E', which is well known as the early integrated Sachs- Wolfe 
(ISW) effect. In the AD and CI modes, the early ISW effect is significant, since their 
gravitational potentials are relatively large at the time of horizon crossing. Contrastively, 
the gravitational potentials are small in the NID mode, and hence their decay affects the 
photon perturbations little. 

3.4 CMB temperature power spectrum 

In Fig. 3, we plot the CMB temperature power spectrum {Cj'^) for the NID mode as well 
as those for other AD and CI modes for reference. The spectrum of the NID mode has 
distinct features coming from the fact that the evolution of photon perturbations in the 
NID mode is different from those in the AD and CI modes as we have seen above. 

First of all, the positions of acoustic peaks in Cj^ for the NID mode are more similar 
to those for the AD mode than those for the CI mode, which reflects the phase of acoustic 
oscillation in the each initial mode. In addition, the amplitude at the first acoustic peak 
(£ ~ 300) in the NID mode is almost same as that at the Sachs- Wolfe (SW) plateau 
(£ < 50). This refiects the smallness of the gravitational potentials in the NID mode, 
where decay of the potentials does not boost the photon temperature significantly before 
and after the recombination. This should be contrasted to the AD mode, where the 
amplitude of Cj^ at the acoustic regime is higher than that at the SW plateau. 

Finally, we would like to give some brief comments on the effects of AA^efi > on the 
CMB power spectrum. At the sensitivity of current and near-future observations such as 
Planck, two characteristic changes from the case of AA'^eg = would be important. One 
is that the RD epoch lasts longer and hence the expansion rate of the Universe becomes 
larger given at fixed redshifts. This enhances the early ISW effect and causes shifts in 
the acoustic peaks towards smaller angular scales. Another is that the free-streaming of 
neutrinos affects the perturbation evolution more significantly. This enhances the decay 
of the gravitational potentials before perturbations enter the sound horizon of the photon- 
baryon fiuid and dampens the photon perturbations in the RD epoch. Thus, this effect 
suppresses the amplitude of Cj"^ at small angular scales. For more detailed discussion, 
we refer to e.g. Refs. [49, 9]. 
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Figure 3: The CMB temperature power spectrum {Cj ) for the AD, CI and NID modes. 

4 Constraints from current observations 

In this section we derive constraints on the isocurvature perturbations in the extra ra- 
diation using current cosmological observations. We adopt CMB data of WMAP 7-year 
result [50, 51, 52] and ACT [53, 54, 2] at small angular scales. While CMB power spectrum 
is the best probe for neutrino isocurvature perturbations, CMB itself contains considerable 
amount of parameter degeneracies among S'dr, N^g and other cosmological parameters. 
In order to solve the degeneracies, we also include data from the baryon acoustic oscil- 
lation (BAO) in the power spectrum of SDSS galaxies [55] and the direct measurement 
of the Hubble constant (HO) [56] . Hereafter, we will refer to sets of combined datasets of 
WMAP+ACT and WMAP+ACT+BAO+HO as "CMB" and "ALL", respectively. 

The initial condition for structure formation in our model is a mixture of the AD and 
NID modes and these two initial modes can be in general correlated. Correlation functions 
of the initial modes, or equivalently the primordial perturbation spectra, form a symmetric 
matrix. We assume that the primordial perturbation spectra can be represented by power- 
law with same spectral indices. Then the primordial power spectra can be parametrized 
as follows, which is widely seen in literatures: 






A.. 



ris-l 



I — a 'y^ya{l 

7a/"(1 -") 



a. 



a 



(96) 



where Vab is the (cross-)power spectrum of initial perturbations A and B defined in 
Eq. (14), and 7 denotes the correlation parameter given in (13). 

In a most general case, the parameter space we explore consists of nine primary cosmo- 
logical parameters (cu;,, tUc, 9s, T, iVeff, Us, As, a, 7). Definitions of these parameters and 
top-hat priors we adopt in the parameter estimation are listed in Table 2. In particular, 
we take A'efr > 3.046 since our model assumes that there are always the Standard Model 
neutrinos which are fully thermalized in the early Universe. In order to take account of 
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UJt 


0.005 -> 0.1 


LOc 


0.01 ^0.99 


9s 


0.5^ 10 


T 


0.01^0.8 


NeS 


(3.046 -> 10) 


Us 


0.5^ 1.5 


in[ioi%; 


2.7^4 


a 


0^1 


1 


(-1^1) 



parameters symbols prior ranges 

baryon density parameter 

CDM density parameter 

angular scale of sound horizon 

optical depth of reionization 

effective number of neutrino species 

spectral index of primordial power spectra 

amplitude of primordial power spectra 

fraction of NID mode in primordial perturbations 

correlation between AD and NID modes 

template amplitude of thermal SZ effect Asz — )■ 2 

template amplitude of Poisson distributed sources Ap 0—7-30 

template amplitude of clustered dust A^i — ?■ 15 

Table 2: Model parameters we constrain in our analysis. First nine from the top are the 
primary cosmological parameters, and rest three are the nuisance ones. The last column 
shows the range of top-hat priors for the parameters we adopt in the analysis. Note that 
the prior ranges of A'^eff and 7, which are shown with parentheses, are only adopted when 
they are varied in the analysis (see text for more details). 



foregrounds, we also include three nuisance parameters Asz, Ap, A^i, which measure the 
amplitudes of power spectra from the Sunyaev-Zel'dovich (SZ) effect, point source, and 
clustered dust, respectively. We adopt the same template power spectra for these fore- 
grounds as in the cosmological parameter estimation of Ref . [2] . In particular, templates 
for the SZ effect and clustered dust are based on Ref. [57]. 

There are several works where constraints on the NID mode are investigated [58, 59, 
60, 62, 63]. A'^eflf is however fixed to the standard value 3.04 in these analyses, so that our 
analysis explore a new parameter space which has not been investigated so far. 

Parameter estimation is performed using a modified version of the publicly available 
CosmoMC code [64]. Convergence of a Markov chain Monte Carlo (MCMC) analysis is 
diagnosed by the Gelman- Rubin test R— 1 < 0.1. 

4.1 Fixed 7 

Before varying all the nine primary parameters, we first fix 7 to or ±1 and vary other 
parameters. This allows us to investigate the constraints for the uncorrelated and totally 
(anti-) correlated cases separately. 

Let us see from the constraints for the uncorrelated case (7 = 0). Constraints on 
parameters are summarized in Table. 3. In particular, we obtain constraints Nq^ < 5.3 
and a < 0.11 at 95 %CL from the ALL dataset. We also show 2d constraints in N^s- 
a plane in Fig. 4. We see there is no strong degeneracy between Ags and a, so that 
constraints on these two parameters can be discussed separately. As is often discussed, 
current CMB observations essentially constrain the epoch of matter-radiation equality. 
Consequently, when Nes is freely varied, Uc and Nek are strongly degenerated. This 
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Figure 4: 68 % and 95 % CL constraints in Ngg-a plane from the CMB (red solid) and 
ALL (green dashed) datasets. From top to bottom, shown are the constraints for the 
uncorrected (7 = 0), totally correlated (7 = 1) and totally anti-correlated (7 = —1) 
cases. Note that the scales are not same among three panels. 
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parameters 


CMB 


ALL 


lOOwf, 


9 qc+O.lO 
^•'-''J-0.10 


9 9cr+U.05 
^■^"^-0.05 


i^c 


n iQc+o.oie 


n 1 QQ+0.012 
U.±OO_Q0i2 


0s 


1 nQS+0-003 

i.UOO_Q003 


1 nQ7+0-003 

±.uor_Qoo3 


T 


0.089+°:°°? 


0.085+°:°S6' 


iVeff (95%CL) 


<6.5 


< 5.3 


Us 


1 02+°°^ 


0.987^°:°}^ 


ln[10^%' 


o 1 q+0.05 


q 1 7+0.04 
'-'•-'-'-0.04 


a (95%CL) 


<0.16 


<0.11 



Table 3: Constraints for the uncorrelated case (7 = 0). We present mean values and 68 
%CL intervals for cosmological parameters except for A'eff and a, for which we present 95 
%CL intervals since they are not bounded from below. 



parameters 


CMB 


ALL 


lOOuJb 


9 1 7+0.06 
^•-■-'-0.06 


2.19+°:°^ 


UJc 


r, 1 qq + 0.014 

U.iOO_Qoi4 


0.130+°:°!^ 


0s 


1 n'?o+0-003 

i.UOZ_QQQ3 


1 n'?'?+0-003 

±.UOO_Q003 


T 


0.087+°:°°? 


0.088i°:rr 


iVeff (95%CL) 


<5.4 


< 5.1 


ris 


U->J'-0.02 


n 979+0-015 
U.J/Z_0.oi5 


ln[10i% 


-| f^+0.04 
'-'•-'-"-0.04 


q 1 c,+0.04 
'-'•-'-^-0.04 


a (95%CL) 


< 0.012 


< 0.011 



Table 4: Same table as in Table. 3 but for the totally correlated case (7 = 1). 



degeneracy is partially solved by including BAO and HO, which determine Vim and Hq 
and hence lead to a better constraint on N^,^. On the other hand, BAO and HO data are 
also effective in the determination of w^, which is strongly degenerated with Ug because 
these two affect the relative height of the first and second peaks in a similar way. a is 
mostly degenerated with ujb and n^ because mixture of the NID mode a 7^ affects the 
relative height of acoustic peaks with little effect on the peak positions. Thus, inclusion 
of BAO and HO improves the constraint on a by solving its degeneracy with Wf, and ns.^° 
Let us then turn our attention to the totally correlated (7 = 1) and anti-correlated 
(7 = — 1) cases. Constraints on parameters for these two cases are summarized in Tables 4- 
5 and 2d constraints in A^eff-a plane are also shown in Fig. 4. From the ALL dataset, 
we obtain A'cfr < 5.4 and a < 0.011 for the totally correlated case and N^fi < 5.2 
and a < 0.007 for the totally anti-correlated case, respectively. We first note that the 
constraints on a for the correlated cases are an order of magnitude tighter than those 
for the uncorrelated case. This is because given fixed a, the NID mode affects the CMB 



-'^"Please notice that if we are allowed to include an extra dark radiation, i.e., A'cff > 3 the best fit value 
for the spectral index increases and the Harrison- Zeldovich spectrum with n^ = 1 is consistent with the 
present observational data. This is true even if the isocurvature perturbation does not exist. 
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parameters 


CMB 


ALL 


lOOW;, 


9 qi+0.07 
^•'-'-'--0.07 


9 9cr+U.05 
^■^"^-0.04 


i^c 


n iQc+o.oie 


n lopr+O.Oll 
U.iOO_o.oil 


0s 


1-041+°:°°^ 


1.039+n^^ 


T 


0.089+°:™^ 


0.083+n^^ 


iVeff (95%CL) 


<6.2 


< 5.2 


Us 


1 01+002 


0.987;ni! 


ln[10^%' 


o 1 Q+0.04 


q 1 1 +0.04 
'-'•-'--'--0.04 


a (95%CL) 


< 0.013 


< 0.007 



Table 5: Same table as in Table. 3 but for the anti-correlated case (7 



power spectrum more significantly for the correlated cases than the uncorrelated one due 
to the power P^^SonW arising from the cross correlation of the AD and NID modes. For 
small enough a of order e.g. 10~^, changes in the CMB power spectrum are somewhat 
similar in magnitude but opposite in the sign. However, data adopted in the analysis 
are not very constraining and relatively large a ~ 0.01 is allowed, where degeneracies of 
a with other parameters are different between the totally correlated and anti-correlated 
cases. For the totally correlated case (7 = 1) a is not degenerated with other parameters 
strongly. This is the reason why the CMB constraint on a is little improved by including 
BAO and HO for this case. On the other hand, a considerably is degenerated with Uc and 
Us for the case of anti-correlated case (7 = — 1). Better determination of Ho and Qm due 
to inclusion of BAO and HO directly tightens the constraint on Uc- Furthermore, n^ itself 
is also strongly degenerated with cUc and hence becomes more tightly constrained. Such 
the improvement in determination of Uc and Ug finally results in better constraints on a. 

4.2 Fixed N^s 

As we have discussed in Introduction, there are several cosmological observations which 
independently show preference for a somewhat larger N^.^ ~ 4 over the standard iVefj = 
3.04 at around 2a level. When we look over the constraints obtained so far, it can be seen 
that such preference still survives, though with somewhat reduced significance, even if 
we include the possibility for the extra radiation to have isocurvature fluctuations which 
can be uncorrelated or totally (anti-)correlated with (. This motivate us to investigate 
a generally correlated case with a fixed N^s 7^ 3.04, for which we devote the rest of this 
section. 

We perform MCMC analysis with 7 being varied and N^s = 4 being fixed. We impose 
a top-hat prior on 7 with a range [—1, 1], and priors for other parameters are as listed 
in Table 2. This time we use only the ALL dataset and the resultant constraints are 
summarized in Table 6. In addition, the 2d constraint in the 7-a plane is shown in Fig. 5. 
As can be seen from the figure, we do not find any evidence for a > and the extra 
radiation which is assumed to exist is consistent with purely AD perturbations. When 
7 is marginalized, we obtain an upper bound a < 0.12 at 95% CL, which is almost the 
same as one obtained for the uncorrelated case. We note that the constraint on 7 is of 
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parameters 


ALL 


lOOwfc 


2 2A+^-^^ 


i^c 


n 1 oq+0.004 


0s 


1 (T^T+O-OOS 

i.uor _o.oo4 


T 


0.086+°:°°^ 


ris 


0.983+°:°}^ 


ln[10i°A; 


3.16+°;°^ 


a (95%CL) 


<0.12 


7 


n 02^+°-^^^ 



Table 6: Constraints in the case that Nes is set to 4 and 7 is varied. We present mean 
values and 68 %CL intervals for cosmological parameters except for a, for which we 
present 95 %CL intervals since it is not bounded from below. 



a 0.1 




Figure 5: 68 % and 95 % CL constraints in 7-0; plane from ALL datasets. 



little importance, as long as a is consistent with zero. 

4.3 Varying N^q and 7 

We finally present constraints for the most general case, where all nine primary parameters 
including Nes, ol and 7 are varied simultaneously. Parameter constraints from ALL dataset 
are presented in Table 7. Compared with the case with fixed A'^gg = 4 in Table 6, we can 
find that there is no significant degradation of constraints on parameters except for cUc, 
which is significantly degenerated with A'eff- On the other hand, upper bound on A"cff is 
also as tight as those obtained for the case of fixed 7 in Section 4.L 

Marginalized constraints on Aeff and a are given as Aeg < 5.1 and a < 0.11 at 95 
%CL. Following the results so far, we again do not find any signature for existence of 
non- vanishing isocurvature perturbations in extra radiations. 2d marginalized constraints 
for A'eff, tt and 7 are presented in Fig. 6. 
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parameters 


ALL 


lOOwfc 


9 9/1+0.05 
^■^^-0.05 


UJc 


r, 101 +0.011 
U.iOi_o.012 


Os 


1 (T^T+OOOS 
±.UO(_ooo3 


T 


0.083+°:°°^ 


iVeff (95%CL) 


< 5.1 


Us 


0.983+°:°}^ 


ln[10^°v4; 


3.1^°^ 


a (95%CL) 


<0.11 


7 


0.07+°:}^ 



Table 7: Constraints for the generally correlated case (—1 < 7 < 1) from ALL dataset. 
We present mean values and 68 %CL intervals for cosmological parameters except for Nes 
and a, for which we present 95 %CL intervals since they are not bounded from below. 



5 Implications 

We have formulated a method for calculating the isocurvature perturbations in the extra 
radiation and derived observational constraints on them. Now let us apply the results 
to some concrete cases. Before discussing a concrete model, we derive formulae for some 
limiting cases for later convenience. 

5.1 Limiting cases 

5.1.1 Uncorrelated case 

First, we suppose that the infiaton does not decay into X {r^ = 1) and the adiabatic 
perturbation is dominantly produced by the infiaton (A^,^ ^ N^-). The a decays into 
X with branching fraction of 1 — r„. Then one can easily imagine that the X has an 
isocurvature perturbation that traces the primordial perturbation of a. 
In this case, from Eq. (58) we have 



3(1 -r^)i?, fg,iH = T,) 



c,(i-(i-r.)i?.) V^*(^ = r 



1/3 



(97) 



and 



with 



^^^ ^ (l-r.)(l-i?.)(l-c.) r g.{H = T,) \ '^^ 2R6a 
Rbr \9*iH = T^)J ai 

RBR^{l-r.)Ra(^^l^^] +c,{l-R„ + r^R„). (99) 

In deriving (98), we have approximated parameters as Rx — Rx ^ 1 and Rr — Rr- Note 
that the condition (A",^ ^ Ao-) leads to the constraint 2R5a /a^ <^ 10~^. In that limit the 
correlation parameter (61) is small : I7I ^ 1. Thus the isocurvature perturbation is nearly 
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Figure 6: 68 % and 95 % CL constraints on (A^eff,a,7) for generally correlated case from 
ALL dataset. 



uncorrelated with the curvature perturbation. It should be noticed that S'dr vanishes in 
the limit Ra ^ 1, since the inflaton contribution to the curvature perturbation becomes 
zero, as is clear from (28). 

Fig. 7 shows constraints on Ra-r^ plane for uncorrelated case. This is drawn by 
translating the results shown in the top panel of Fig. 4 into the constraint on R„ and r^, 
using (58) and (59). la and 2a allowed regions are shown by blue and orange, respectively. 
Contours of a and ANes are also shown. In this figure 2R6a/ai = 0.1x(^^~5x 10~^ 
and g*{H = T„) = 100 are fixed. 

5.1.2 Totally anti-correlated case 

Next, suppose that the a does not decay into X, i.e., r„ = 1, and a dominantly produces 
the curvature perturbation (A^^ ^ A^^-) : o" truly takes a role of curvaton. The inflaton 
is assumed to decay into X with branching fraction of 1 —r^. Then the extra radiation X 
produced by the inflaton decay is expected to have isocurvature perturbation correlated 
with the curvature one. 



29 



0.8 



0.6 



0.4 



0.2 




Figure 7: Constraints on Ra-Ta plane for uncorrelated case. Icr and 2a allowed regions 
are shown by blue and orange, respectively. Contours of a and AA^eff are also shown. In 
this figure 2R5a/ai ~ 5 x 10~^ is fixed. 



In this case we have 



AA/", 



cff 



3(l-r^)(l-i?,) [g,{H = T, 



Cu{r^{l - R„) + R„) \g^{H = T^ 



1/3 



and 



with 



5i 



DR 



(l-r^)(l-i?.)(l-c,) fg.iH = T 



R 



DR 



g*{H = T^) 



1/3 



2R6a 



CTi 



Rdr ^ (1 - r^){l - Ra) [ ^*^Z ^i 1 + c-(^<^ + ^- - ^^^-) 



(100) 

(101) 
(102) 



In deriving (101), we have approximated parameters as Rx — Rx ^ 1 and Rr ~ Rr- 
Notice that we have the condition 2R6a/ai ~ 1.5 x 10~^ in order for the curvaton to 
produce the curvature perturbation. In this limit, the correlation parameter (61) becomes 
7 ~ — 1 : the isocurvature perturbation is almost totally anti-correlated with the curvature 
perturbation. Also we need R > 0.01 in order not to have too large non-Gaussianity. 

Fig. 8 shows constraints on R^-t^ plane for totally anti-correlated case. This figure is 
obtained by converting bottom panel of Fig. 4 using (58) and (59). 2a allowed regions 
are shown by orange. Any values of R^ and r^ are excluded at la level. In this figure, 
contours of both a and AA'eff are shown. Note that in the case of To- = 1, the following 
relation between a and AA^gg is satisfied. 



a 



S^ 



DR 



3(1 - c^ 



AK 



cff 



(103) 



1-a ( "'" ""'S + AA/'eff 

which is derived from (58), (59) and (60). Here g*{H = Fg-) is set to 100 in this figure. 
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Figure 8: Constraints on R^-r^ plane for totally anti-correlated case. 2a allowed regions 
are shown by orange. There is no la allowed region. 



5.2 Particle physics model 

A model to explain AN^s ~ 1 was proposed in Ref. [17]. It is based on the supersymmet- 
ric (SUSY) extension of the axion model. The axion is a pseudo Nambu-Goldstone boson 
associated with the spontaneous breakdown of the Peccei-Quinn (PQ) symmetry, which is 
introduced in order to solve the strong CP problem in quantum chromo dynamics [65] . In 
a SUSY axion model, there exists a scalar field a, called saxion, which is a scalar partner 
of the PQ axion [66] . The saxion can naturally have a mass of the order of the gravitino, 
which is the superpartner of the graviton, and it can be 0{keV)-0{TeV). The saxion can 
have a large initial amplitude ai during inflation, and may obtain quantum fluctuations 
of 6a ~ ifinf/(27r) since its mass is much smaller than the Hubble scale during inflation. ^^ 
Therefore, the saxion is a candidate for the curvaton. The saxion decay mode depends 
on the model (see e.g. Refs. [67, 68, 69, 70, 71, 72]). We consider following two cases : 
the saxion dominantly decays into two axions {a — t- aa) and into the Higgs boson pair 
{a —7- hh). The former situation is typically realized in the KSVZ axion model [73] and 
leads to the uncorrelated isocurvature perturbation in the relativistic axion particles if 
the curvature perturbation is dominantly sourced by the inflaton. The latter possibility 
appears in the DFSZ axion model [74]. In this case, the saxion can take a role of the 
curvaton and thermally produced axions after the inflaton decay have correlated isocur- 
vature perturbations. In the axion model, the axion also has an isocurvature perturbation 
which results in the CDM isocurvature mode. We do not discuss it in detail since it can 
be reduced by tuning the initial misalignment angle. 



^^The Hubble-induced mass term of the saxion is assmiied to be suppressed. Otherwise, its quantum 
fluctuation is significantly suppressed. 
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5.2.1 SUSY KSVZ axion model 

First, let us consider the case where the saxion decays into two axions. This is often the 
case in the KSVZ axion model [73] . The decay rate of this process is given by 

where fa denotes the PQ symmetry breaking scale, which is constrained as lO^GeV ^ fa '^ 
lO^^GeV and / is a typically 0{1) constant. Hereafter we take / = 1. The temperature 
at the saxion decay is calculated as 

/ 100 Y^^ f m, \3/2 /ioi2GeV\ , , 

T„ - 100 MeV — — - ( -f— ) . 105 

V^*(^-)/ VlOOGeV/ \ fa ) ^ ' 

The saxion also decays into gluons if kinematically allowed. The branching ratio of the 
saxion into gluons, r^, is given by 



2a 



r^ ~ — f ~ 3 X 10"^ (106) 

and hence very small. Nonthermal axions produced by the saxion decay is regarded as an 
extra radiation X, since the PQ axion has only extremely weak interaction with ordinary 
matter and its mass is very small. The abundance of the saxion coherent oscillation is 
given by 

'^ VKf^V^2xlO-Gevf,,^V,,,,STT.,)^ (107) 



s 8 "\MpJ VlO^GeVy Vl0i2GeV 

for rria > T^ and 

'- = -T.3C f ^ V - 2 X 10-^GeV ( "^^ ) ''' ( ,"- ) ' , (108) 

s 8 °''V^P/ VlOOGeV/ Vl0i2GeV/ ' ^ 

for rricr < T^ where s denotes the entropy density, Tr the reheating temperature after 
inflation which is related to the infiaton decay rate P^ as Tr = (10 / n'^ g^Y^^ ^yT ^Mp and 
^osc = (10 /Tr'^g^Y^^y/rriaMp. Thus R ~ (3/4)_Ro- = 3pcr/(4ptot) at the a decay is estimated 
as 

U=l^^^2.10-'(J^](l^r(,J^X('-^)\ (109) 



for i? ^ 1 and nia- > P,/, and 

i? - i^^ ~ 2 X 10-^ (^-^^] ( ^" ^ ' (^X (110) 

for i? ^ 1 and ma < P</,. Assuming that the infiaton decays only to the visible sector 
{r^ = 1), ANes and Sbr are obtained by substituting (109) into Eqs. (97) and (98). The 
AXeflF is estimated as 

Therefore, both ANcs and Sdr can have sizable values depending on model parameters. 



32 



5.2.2 SUSY DFSZ axion model 

Next let us consider the case where the saxion doniinantly decays into a Higgs boson pair, 
as is realized in the DFSZ axion model [74]. In this case, the saxion can take a role of the 
curvaton. The decay rate of this process is given by 

where /i denotes the so-called /x-parameter in SUSY, which gives the higgsino mass. If 
this is the dominant mode, the invisible branching ratio into an axion pair is estimated 
as 

1 /m„ 



which can be very small if /i > ma, and we have r^ — 1. However, even if we neglect 
the nonthermal axions produced by the saxion decay, there is another contribution from 
thermal production during reheating. For simplicity we consider the case that axions are 
thermalized after the infiaton decays, which occurs if the reheating temperature satis- 
fies [75] 

/ f \ 2-246 

T^>T^^ 9.6 X 10^ GeV [^0^) ' ^^^^^ 

If the reheating temperature is higher than this, axions are thermalized for T > Td and 
decouple from thermal bath at T ~ Td. At this moment, the fraction of the axion energy 
density relative to the total radiation energy density is given by 5'*(T = Tq)^^. Thus 
we can effectively regard that the infiaton decayed into axions with branching fraction 
of ~ 5f^,(T = Td)~^. Therefore, we can replace 1 — r,^, which is defined at the curvaton 
(saxion) decay, with 

From Eq. (100), we obtain a well known expression for the AA'^gg for thermalized species 
(see e.g., Ref. [19]) 



g^{H = Ta)J 7 \g^{T = Ty) 

for Ra <^ 1. Otherwise, the saxion decay dilutes the thermally produced axions to a 
negligible level. Those thermal axions have an anti-correlated isocurvature perturbation 
to the curvature perturbation, since the latter is assumed to be produced dominantly by 
the curvaton. The magnitude of the isocurvature perturbation in this case is given by 
(101). 

6 Conclusions and discussion 

In this paper, we have investigated isocurvature perturbations in an extra radiation com- 
ponent. We have firstly formulated how primordial isocurvature perturbations between 
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the extra radiation and the Standard Model radiation can be generated from fluctua- 
tions of two scalar fields in the inflationary Universe. We gave a detailed description 
of what isocurvature perturbations in the dark radiation are generated. Our derivation 
of the isocuravature perturbations is based on the SN formalism and fully non-linear. 
We have also pointed out that non-Gaussianities can be generated in such isocurvature 
perturbations, which are to be investigated in more detail in future works. It would be 
straightforward to extend our formulation for cases with three or more fields. We have 
also discussed observational signatures of the isocurvature perturbations in the extra ra- 
diation. We have pointed out that our model leads to distinct features in the CMB 
power spectrum which are different from those predicted by the ordinary curvature and 
CDM/baryon isocurvature perturbations. We have shown that that our model can be 
constrained from current cosmological observations. Roughly speaking, CMB combined 
with BAO and the direct measurement of Hq gives constraints on the abundance of the 
extra radiation and the amplitude of its isocurvature perturbations as 3 < iVeff < 5 and 
a < 0.1 (0.01) for the uncorrelated (totally correlated/anti-correlated) case. These are 
also translated into constraints on some limiting scenarios which can be realized in models 
based on particle physics such as SUSY axion model. 

Our model will be tested more precisely by cosmological observations which are on- 
going or projected in the near future. In particular, CMB power spectrum from the 
Planck survey is expected to improve constraints on both N^s and a by an order of mag- 
nitude. Furthermore, non-Gaussianities and other observational signatures would be also 
informative. 
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A Extra radiation production after BBN 

So far we have assumed that the a decays before the BBN and hence the neutrino freezeout 
epoch. This may not necessarily be the case. If the a mainly decays into the X particle 
without producing a significant amount of visible particles, the decay of a does not upset 
BBN. Thus it may be possible that the a decays after the neutrino freezeout but well before 
the recombination. In this case the produced X particles contribute to AA'^eg measured 
from CMB but does not affect ANes measured from BBN. Examples of concrete models 
were proposed in Ref. [17, 20]. In this appendix we derive a formula in such a case. 

The scenario we are considering here is summarized in Table 8. The curvaton a decays 
after BBN begins, if the decay is too late, the effects of extra radiation on the CMB and 
LSS are much more different and we do not consider such a case. Thus we restrict the 
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lifetime of the curvaton to be less than ~ 10^ sec. By taking the uniform density slice at 
the curvaton decay, we have following relations 

pf^S) + (1 _ ,(7) _ r('^))p^(f) = p^(f), (117) 

pf(f) + r(V(^)=P7(^), (118) 

pW(f)+r(^W(^)=P.(:?), (119) 

pf{x) + pW(x) + pf{x) + p.(f) = ptot(x)(= ptot). (120) 

Here r^ and r^ are the branching ratios of a into the photon and neutrino. ^^ The 
curvature perturbation of each component is related to its background value as 

(121) 

(122) 



pf(^) 



^We4(c,-C)^ 



pW(f) 






p-e-- -, (123) 

Pa{x) = p^e^K"--^). (124) 

The requirement that the curvaton decay should not spoil BBN sets strong constraint 
on the branching ratio into the photon [76] and neutrino [77] . Therefore we can safely set 
Ta = Tct = since we are considering the case where the curvaton energy density is not 
negligible at the curvaton decay. Then we soon find 

C7 = a, Cu = C^. (125) 

The total curvature perturbation and the curvature perturbation in X are found to be 



+ i?(C. - C0) + 2^(1 - ^)(3 + Rm - Ca) 



c 



X 



m 



X 



?(0) 



5{^) 



(47?w + (1 _ R)B}^>)C^ + (3 + R)B}^Xa 



W. 



(126) 



+ 



(3 + R)Rf^ 



"-x 



^R-x 



M 



(^) 



l + R){3-R)R'^> + {l-R)RR'^> (^ - C<x) 



?(^) 



where R is defined by Eq. (29) with R^ = Pa/ptot\H=ra- Notice also that R^' = R^ and 
Rr = in the present case. Other quantities are also defined as Eq. (27). Since there are 
no changes in the relativistic degrees of freedom after the curvaton decay, the expression 
(126) gives the final curvature perturbation. The DR energy density is defined by 

p(^)(x) + pf{x) + p^(f) = pdr(x), (127) 

on the uniform density slice, where Py)k{^) = Pdrc'^^''^'^"''-*. The DR isocurvature pertur- 
bation is calculated as 



^1 



DR 



3(Cdr - C) 

-3^(1- 
-Rdr 



-(C-C.) + ^ ^^-^f^^^ £(l + i^OR) 



{Co 



fl28l 



^^The decay products of a, except for the neutrino and X, are thcrmalized. They are coUectively caUed 
"photon" here. 
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epoch 


component 


energy transfer 


T4,<H 


0, a 


^ XW + rW 


V,<H <V^ 


XW, rW, a 


j^W _^ ;y + ^^ 


T^± <H < Ty 


XW, z/, re, a 


e=^ -)■ 7 


T,<H < Tei 


XW, z/, 7, a 


a -^ X^'^^ 


H<T, 


X, I/, 7 (DR = X + z/) 





Table 8: Same as Table 1 but for the case of T^ < Fgi. The a is assumed to dominantly 
decay into X. 



where Cj. is given in Eq. (51), Rr = (p^ + Pu)/ptot and Rdr = {px + Py)/ptot evaluated 
after the curvaton decay. 

The extra effective number of neutrino species, AXefj, is given by 



AXeff = 3 

This is exphcitly evaluated as 

AXeff ^ 1.1 



Px _ 3-R<7 

Pu Ci/lL'p 

/lkeV\ / Pa/s 



(129) 



/lkeV\ / pjs \ 
\T, ) VlO-^GeVy' 



(130) 
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